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Abstract 

Many  modern  control  design  methods,  such  as  H2  and  H . 
optimization,  guarantee  that  the  closed  loop  system  will  be 
stable.  To  make  further  use  of  this  guarantee,  the  portion 
of  the  s-plane  considered  to  be  stable  is  restricted  to  a 
region  left  of  a  vertical  line  some  distance  a  from  the 
imaginary  axis.  This  a  shift  was  accomplished  in  the  state 
space  equations,  so  standard  computational  tools  still  work. 
Now  the  system  is  not  only  guaranteed  to  be  stable,  but  also 
has  some  guaranteed  speed  for  the  time  response.  The 
equations  for  guaranteeing  some  arbitrary  damping  ratio  were 
also  derived  in  Cartesian  coordinates,  but  only  the  a  shift 
was  put  into  the  state  space  equation  form  and  tested.  The 
a  shift  increased  the  speed  of  the  system,  as  expected.  An 
additional  benefit  was  increasing  the  damping  of  lightly 
damped  poles.  The  cost  for  doing  this  was  an  increase  in 
the  steady  state  error  and  an  increase  in  the  control  power 
required.  The  a  shift  is  very  easy  to  perform  in  the  state 
space  equations.  The  bottom  line  is  that  it  gives  the 
design  engineer  another  parameter  to  adjust  to  get  better 
system  performance. 
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H 2  AND  H m  OPTIMIZATION  WITH  A 
RESTRICTED  REGION  OF  STABILITY 


1 >  Introduction 


Many  modern  control  design  methods  guarantee  stability, 
i.e.  that  all  the  closed-loop  poles  are  in  the  Left  Half 
Plane  (LHP).  Since  stability  is  insured  by  the  mathematics 
of  the  design  procedure,  the  design  engineer  can  focus  his 
attention  on  achieving  other  design  parameters.  Rather  than 
just  guaranteeing  stability  it  would  be  beneficial  to  insure 
other  parameters,  such  as  a  minimum  speed  of  response  and  a 
minimum  damping  ratio.  This  could  be  accomplished  by 
restricting  the  region  of  the  s-plane  considered  to  be 
stable.  This  thesis  will  present  a  method  for  transforming 
the  "stable"  region  of  the  s-plane  to  an  area  left  of  a 
vertical  line  some  distance  a  to  the  left  of  the  imaginary 
axis.  Shifting  all  the  closed-loop  poles  away  from  the 
imaginary  axis  will  increase  the  speed  of  the  system  and 
will  increase  the  damping  of  complex  poles  that  are  close  to 
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the  imaginary  axis. 

Chapter  2  presents  the  equations  and  design  methodology 
for  H2  and  H .  optimization.  It  also  gives  some  of  the 
properties  and  uses  for  each.  Chapter  3  gives  an  overview 
of  the  previous  work  in  this  area.  Since  H2  and  Hm 
optimization  are  so  recent,  the  earlier  work  was  for 
designing  Linear  Quadratic  Regulators  (LQR) .  Chapter  4 
presents  the  mathematics  for  transforming  the  complex  plane. 
Besides  the  a  shift,  it  shows  the  equations  for  rotating  the 
imaginary  axis  by  some  angle  This  transformation  will 
guarantee  some  minimum  damping  ratio.  Chapter  5  shows  the 
state  space  representation  of  the  a  shift.  Chapters  6  and  7 
are  example  problems.  The  first  one  is  a  simple  plant  that 
shows  the  effect  of  shifting  the  imaginary  axis  across 
open-loop  poles  and  zeros.  The  second  example  is  for  a 
drone  aircraft.  It  shows  how  the  a  shift  can  be  used  to 
increase  the  damping  of  a  system  with  complex  poles  near  the 
imaginary  axis.  Chapter  8  shows  how  the  a  shift  affects 
robustness.  The  compensator  that  was  designed  for  the 
original  plant  is  tested  on  several  perturbed  plants. 

Chapter  9  gives  the  conclusions  and  recommendations  for 
areas  of  further  study. 
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2 .  H 2  and  H .  Control  Design  Methods 


Classical  control  design  methods,  such  as  root  locus  ard 
Bode  plots,  work  very  well  on  simple  Single  Input  Single 
Output  (SISO)  systems.  As  the  system  becomes  more  complex, 
however,  new  design  methods  are  needed.  For  a  complex  SISO 
system  simply  feeding  back  the  gain  or  adding  filters  may 
not  be  enough  to  achieve  the  desired  closed-loop  response 
and  for  Multiple  Input  Multiple  Output  (MIMO)  systems  the 
classical  methods  are  not  easy  to  use  at  all.  The  following 
discussion  on  modern  control  design  methods  is  not  meant  to 
be  complete  or  quantitative.  The  interested  reader  is 
referred  to  [4],  [5],  [8]  and  [10]  for  complete  discussion 
and  mathematical  derivations.  H2  and  Hm  optimization  will 
be  the  only  methods  discussed  here.  Linear  Quadratic 
Gaussian  (LQG)  is  in  the  category  of  modern  control ,  and  is 
a  special  case  of  Hz  optimization. 

2.1.  B1 ock  Diagram  and  State  Space  Equations 
Consider  the  block  diagram  in  Figure  1. 
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Figure  1  General  feedback  block  diagram 

The  state  space  equations  for  the  system  are 

x  -  Ax  +  Bdd  +  Buu  (1) 

e-  C.x+  D,dd+  Dtuu  (2) 

y  -  Cyx+ Dydd+ Dyuu  (3) 

The  following  are  also  required; 

(A,Bd)  stabilizable  and  (Ct, A)  detectable  (4) 

stabilizable  and  (C,,.d)  detectable  (5) 
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Also,  assume  that 


D«i  “  0  and  Dyu  m  0 
Dr.uDtu-  I  and  DydDTyd-l 


(6) 

(7) 


To  satisfy  the  assumption  in  Equation  (7), 

be  scaled.  The  rescaled  equations  are 

y  and  u  must 

u-s;‘a 

(8) 

y  -  Syy 

(9) 

where  Su  and  Sy  are  found  from 

■SttSu  ■  £>.uD.u 

(10) 

~  DydDTyd 

(ID 

by  using  a  Cholesky  decomposition.  The  scaled  state  space 
equations  become 
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x  -  Ax  +  Bdd  +  Bua 


(12) 


e-Ctx+D.aCl  (13) 

?-Cyx  +  Dydd  (14) 


where 


fl.'fi.s:1  Cy-SyCy  (15) 

fi»u.~DtuSu  DydmSyDyd  (16) 

The  above  notation  will  be  used  for  the  remainder  of 
this  chapter.  What  follows  is  known  as  the  Doyle/Glover 
parameterization  of  H2  and  Hm  optimal  compensators. 

2 . 2 .  Properties  and  Equations  of.  H z  and  H .  Optimization 
2.2.1.  Ho  Optimization. 

H2  optimization  begins  with  the  family  of  all 

stabilizing  compensators  K,  and  finds  the  one  that  minimizes 
the  2  norm  of  the  closed-loop  transfer  function,  Ttd.  The  2 
norm  is  defined  as 
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Ct-v-D'.Dlyc, 


(23) 


and 


(A- BdDTydCy)Y2+Y2(A- BdDTydCy)T -Y2CTyCyY2+  BXB\-Q 


(24) 


where 


Bl-Bd{l-DT„t>^'>  (25) 

Another  way  to  think  of  H2  optimization  is  to  minimize 

the  area  under  the  Singular  Value  (SV)  plot  of  the 
closed-loop  transfer  function.  Obviously,  the  SV  plot  must 
approach  0  as  caj  approaches  »  or  the  2  norm  will  be 
infinite.  This  is  the  reason  that  the  D,d  must  be  0 .  If  it 
is  anything  but  0  the  closed-loop  system  will  not  be 
strictly  proper,  and  the  SV  plot  will  not  approach  0  as  the 
frequency  approaches  <».  The  H2  control  design  method  is 
most  effective  on  a  system  with  a  white  noise  input.  It 
will  design  a  compensator  K  that  minimizes  the  energy  of  the 
output  based  on  a  white  noise  input  and  the  amount  of 
control  power  allowed.  H2  optimization  is  also  good  for 
reducing  the  effects  of  parameter  variations. 
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2.2.2.  H -  Optimization. 


//.  optimization  begins  with  the  family  of  all 

stabilizing  compensators  K,  and  finds  one  that  minimizes  the 
00  norm  of  the  closed-loop  transfer  function,  7^.  Unlike  Hz 
optimization,  the  solution  to  Hm  optimization  is  not 
unique.  The  «>  norm  is  defined  as 

\T*\.m  sup  ommx[T,d(ju)')]  (26) 

tjj 

where  0max  is  the  maximum  singular  value. 

The  central  member  of  the  family  of  all  stabilizing 
compensators  such  that 


(27) 


is  given  by 


K  ~  -  K c(s I  -  A kYl  K  f 


(28) 


where 
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Ak  -  A  -  KfCy  -  BUKC  +  y-2Y.cUct  -  DtuKc ) 


(29) 


(30) 

Kf-y.Cl  +  BaDT,t 

(31) 

X .  and  Ym  are  the  solutions  to  the  following  Riccati 
equations : 


(y4-£u£Lor*-+*-M-Btt/5lc.)-x.(8u£i-v‘2£dSd:)X-  +  cic,«o 


(32) 


where 


c,  -(/-/>  .u^L)C. 


(33) 


and 


(34) 


where 


(35) 
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For  this  to  be  a  valid  parameterization,  the  following  three 
conditions  must  be  met: 


X 

1 

rv 

o 

(36) 

y.2t  o 

(37) 

p  (r-X.)^v2 

(38) 

where  p(-)  is  the  spectral  radius. 

The  Hm  optimization  equations  must  be  iterated.  If  any 

of  the  above  conditions  fail,  y  must  be  increased.  If  all 
of  them  pass  y  can  be  decreased.  By  doing  this  repeatedly, 

Y  can  be  brought  arbitrarily  close  to  the  optimal  value. 

From  the  definition  (Equation  26)  it  is  clear  that  the  « 

norm  is  simply  the  peak  value  of  the  maximum  SV  of  the 
closed-loop  transfer  function.  optimization  has  two 

main  uses.  First,  it  will  minimize  the  energy  of  the  output 
to  an  unknown,  but  bounded  energy  input.  Second,  it  will 
minimize  the  effects  of  unmodelled  dynamics  and  structural 
uncertainties,  i.e.  robustness. 
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3.  Previous  S-plane  Transformations 


The  previous  chapter  discussed  some  of  the  advantages  of 
the  H2  and  //.  optimization  design  methods.  One  of  the 
important  properties  is  their  guaranteed  stability.  Many 
researchers  have  tried  to  expand  on  this  point.  If  all  the 
closed-loop  poles  can  be  forced  to  be  left  of  the  imaginary 
axis,  why  can't  they  be  restricted  to  some  subset  of  the 
left  half  plane?  This  chapter  will  examine  some  of  the 
schemes  put  forward  to  make  further  restrictions  on  the 
location  of  the  poles.  Basically,  the  restrictions  are  one 
of  two  kinds.  First,  move  all  the  poles  some  minimum 
distance  from  the  imaginary  axis;  i.e.  left  of  a  vertical 
line.  Second,  move  only  the  lightly  damped  poles  away  from 
the  imaginary  axis,  or  at  least  move  them  more  than  the 
poles  with  better  damping.  The  first  kind  will  increase  the 
stability  and  speed  of  the  system  and  the  second  kind  will 
increase  the  damping  ratio.  Previous  works  have  examined 
the  theory  for  restricting  the  region  of  stability  for  the 
Linear  Quadratic  Regulator  (LQR)  design.  The  following  is  a 
brief  description  of  some  of  their  techniques. 
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3.1.  Closed-Loop  Poles  Left  of  a.  Vertical  Line 


The  first  effort  to  restrict  the  region  considered  to  be 
stable  was  put  forward  by  Anderson  and  Moore  in  1969  ([2] 
and  [3]).  To  make  it  easier  to  follow,  their  equations  have 
been  put  into  standard  notation. 

They  begin  with  the  time-invariant  linear  system 

X  -  Ax*  Bu  (39) 

u--Krx  <4°> 

subject  to  the  following  performance  index 


J -  f  (urRu+ xtQx 

•'to 


)cit 


(41) 


where  R  is  a  positive  definite  symmetric  constant  matrix  and 
Q  is  a  nonnegative  definite  symmetric  constant  matrix. 

To  get  the  desired  pole  movement  they  replace  the 
performance  index  with 


J -  f  (uT/tu  +  xTQx)e2a‘dt 
J  «o 


(42) 
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where  a  is  the  desired  minimum  distance  left  of  the 
imaginary  axis.  Next  let 


j?  -  xeat 

i2  -  uea‘ 

Substituting  into  the  performance  index 


(43) 

(44) 


J 


(  CLrRa+  xTQ£')dt 


(45) 


If  the  eigenvalue  corresponding  to  any  state  x  does  not 
decay  at  least  as  fast  as  e~at,  J  will  be  infinite.  In  other 
words,  the  real  part  of  each  closed- loop  eigenvalue  must  be 
left  of  -a. 

The  new  state  equation  becomes 

+  (46) 
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The  new  state  equation  is  very  similar  to  the  original. 

A  is  replaced  with  A  +  al  and  everything  else  remains  the 
same.  Anderson  and  Moore  then  examined  some  of  the 
properties  of  the  shifted  system. 

In  1985  Amin  [1]  used  the  principle  discovered  by 
Anderson  and  Moore,  along  with  pole  placement  techniques,  to 
devise  a  method  of  shifting  the  real  parts  of  the  open-loop 
poles  to  any  desired  position  without  affecting  the 
imaginary  parts.  His  approach  moved  each  pole  or  complex 
pair  individually  and  could  be  used  iteratively  to  move  any 
number  of  poles.  His  solution  was  based  on  solving  either  a 
first  or  second  order  Lyapunov  equation;  no  (non-linear) 
algebraic  Riccati  equations  are  involved.  This  approach 
would  be  useful  if  it  is  important  that  the  imaginary  part 
of  the  pole  not  change. 

In  1986,  Sheieh,  Dib,  and  Mcinnis  [12]  presented  a 
method  for  restricting  the  closed-loop  poles  to  a  vertical 
strip  in  the  LHP.  They  began  with  the  standard  system, 
performance  index  and  control  law  given  by 
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x - Ax+ Bu 


(47) 


J  -  f  Qur  Ru  + xrQx)dt 
Jo 


(48) 


u--Kx-~R~lBT  Px 


(49) 


where  P  is  the  solution  to  the  algebraic  Riccati  equation 


PA  +  ATP  +  Q-PBR~lBTP-  0  <5°) 

The  optimal  closed-loop  system  is 

x-(4-5/Ox  (51) 

They  showed  that  the  closed-loop  system  can  be  designed  to 
be 


x  -  {A  -rBK~)x 


(52) 


where 
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(53) 


1  i  h.2~  h  i 
r~2*  tr(BK') 

and  hy  and  h2  are  the  distances  from  the  imaginary  axis  to 

the  left  and  right  sides  of  the  vertical  strip, 
respectively.  This  method  not  only  moves  all  the  poles  left 
of  some  line,  but  also  restricts  them  from  being  too  far  to 
the  left.  It  would  be  useful  on  systems  with  two  time 
scales,  one  fast  and  one  slow. 

3.2.  Closed-Loop  Poles  Inside  &  Cone 

In  1983,  Kawasaki  and  Shimemura  [7]  designed  an 
iterative  method  to  choose  the  quadratic  weights  in  LQR. 

They  approximated  a  45  degree  cone  in  the  LHP  with  the 
hyperbola  (see  Figure  2) 

(/ee\)2-(/m\)2-m2  <54> 

Using  the  standard  state  equation  and  quadratic  cost 
function  they  first  solve  the  standard  LQR  problem. 
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Figure  2  Approximating  a  conic  sector  in  the  LHP  with  a 

hyperbola 

Next,  beginning  with  the  standard  solution,  they  iterate  the 
following  equations  until  Kt  is  equal  to  zero.  Let 

A,m  At_, -BR'xBtP]  <55) 

where  i  =  1,2,...  and  A0m  A 

Next,  solve  the  following  equation  for  £<  (the 
superscript  +  denotes  the  maximum  solution) 
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KlBlTiBTKl  +  Ki('4?-m2I)  +  ('4?-rn2I)TKl-0 


(56) 


If  •  0,  then  the  solution  is  complete.  If  #  0  then 

choose  some  arbitrary  real  number  rt  such  that  rt>  \  an<^ 
calculate  P%\  of  the  following  equation 

Pl.lBR-'BTPl.l-Pl.iAi-ATlPl.x-  rtK;-0  (57> 

Next,  solve  for  Altl  -  At  -  BR~l  BT  P*t+X ,  the  closed-loop  system 

matrix,  and  go  back  to  Equation  (55). 

The  final  closed-loop  system  is 

A- BR~l  Bt(P\+  P*2+  ...+  P])  <58> 

where  P~,  is  the  solution  to  Equation  (57)  for  K\m  0. 

The  above  closed-loop  system  has  all  it's  eigenvalues  in 
the  hatched  region  of  Figure  2. 

Next,  in  1988,  Sheieh,  Dib,  and  Mcinnis  [11]  presented  a 
method  to  restrict  the  closed-loop  poles  to  a  cone  with 
either  a  45  degree  or  30  degree  angle  from  the  negative  real 
axis.  Again  they  begin  with  the  standard  state  equation  and 
quadratic  cost  function.  Their  method  is  an  iterative 
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method  similar  to  the  one  presented  by  Kawasaki  and 
Shimemura.  The  details  and  specific  design  steps  are  in 
[11]. 

One  additional  paper  is  worth  mentioning.  In  1984, 
Oliver  [9]  presented  a  scheme  for  restricting  the  region  of 
stability  to  any  region  symmetric  about  the  real  axis. 

While  this  method  is  general,  it  would  be  hard  to  implement 
in  anything  other  than  very  simple  cases. 
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4.  Plane  Transformations 


4.1.  a  Shift 

Hz  and  H .  optimization  are  extremely  powerful 

compensator  design  methods.  They  guarantee  that  the 
closed-loop  system  is  stable  and  minimizes  either  the  2  norm 
or  the  oo  norm,  respectively,  of  the  closed-loop  transfer 
function.  The  stability  guarantee  is  particularly  useful. 

It  says  that  all  the  closed-loop  poles  will  be  to  the  left 
of  the  imaginary  axis.  For  example,  a  plant  with  a  pole  and 
a  zero  in  the  right  half  plane  (unstable  and  non-minimum 
phase)  will  become  closed-loop  stable.  Regardless  of  where 
the  pole  and  zero  are  in  the  right  half  plane,  Hz  or  Hm 
optimization  designs  a  compensator  that  places  all  the 
closed-loop  poles  in  the  left  half  plane.  Rather  than  just 
placing  all  the  poles  in  the  left  half  plane,  it  may  be 
useful  to  place  them  to  the  left  of  a  vertical  line  some 
distance  a  from  the  imaginary  axis.  This  will  insure  that 
the  closed-loop  system  is  not  only  stable,  but  has  some 
specified  degree  of  stability.  Looking  at  this  from  the 
time  domain,  it  will  insure  that  the  system  has  some 
guaranteed  speed  of  response.  Figure  3  shows  the  a  shift  in 
Cartesian  coordinates.  For  the  equations  in  this  section 
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the  variable  with  the  bar  above  it  is  the  transformed 
variable.  For  example,  x  is  the  original  coordinate  and  x 
is  the  transformed  coordinate. 


Figure  3  Alpha  shift 

The  equations  for  this  plane  transformation  are  simply 


y-y  (59) 


x-x  +  a  (60) 
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The  next  logical  step  is  to  transform  these  equations  to 
the  s-plane  and  finally  to  state  space  equations.  Almost 
all  modern  control  design  calculations,  including  H2  and  Hm 
optimization,  are  performed  on  state  space  equations.  The 
next  chapter  will  show  the  manipulations  required  to 
transform  the  state  space  equations  to  accomplish  an  at 
shift.  Before  that,  however,  some  other  useful 
transformations  on  the  s-plane  stable  region  will  be  shown. 

4.2.  £  Rotation 

As  shown  in  Section  4.1,  transforming  the  left  half 
plane  to  a  sub-plane  left  of  a  vertical  line  some  distance  a 
from  the  imaginary  axis  is  not  very  difficult.  Another 
useful  transformation  would  be  to  rotate  the  imaginary  axis 
by  some  angle  into  the  left  half  plane.  Forcing  all  the 
closed-loop  poles  to  lie  in  this  restricted  region  will 
guarantee  that  the  system  will  have  some  minimum  damping 
ratio. 

4.2.1.  Cylindrical  Coordinates. 

The  first  step  is  to  show  the  rotation  in  cylindrical 
coordinates.  For  this  transformation,  r  (the  distance  from 
the  origin)  is  held  constant.  Figure  4  shows  the  rotation 
of  the  imaginary  axis,  y  is  the  angle  from  the  negative 
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real  axis  to  the  desired  imaginary  axis  location  and  6  is 
the  angle  from  the  positive  real  axis  to  the  point  of 
interest . 


y 


Figure  4  Rotating  the  imaginary  axis  to  some  damping  ratio 


The  following  equation  relates  y  to  the  damping  ratio  £ 


y  "  cos  '1  £ 


0  is  the  angle  from  the  positive  real  axis  before  the 
transformation  and  0  is  the  angle  after  transformation. 
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If  0  =  180°  -  y,  then  0  =  90 


If  0=  0°,  then  0=0° 

If  0  =  180°,  then  0  =  180° 

Figure  5  shows  how  0  is  related  to  0.  Each  line  at 
angle  0  from  the  positive  real  axis  is  rotated  to  some  angle 
0  by  the  following  equations. 


0  £  180° -  y  =>  0- 


90°  „ 

- 9 

180* - y 


(62) 


0>  180°  -  y  -> 


_  90 # 

0-— (0-(18O°-y))^9O# 


(63) 


Simplifying  this  equation  yields 


0  >  180*  -  y  -> 


90° 

0-  180°+ - (0-  180°) 

y 


(64) 
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The  next  step  is  to  convert  these  equations  to  Cartesian 


coordinates.  For  this  transformation,  the  radius  between 
coordinates  is  assumed  constant,  i.e.  rotate  the  damping 
ratio  line  to  the  imaginary  axis  keeping  the  distance  from 
the  origin  constant.  The  equations  to  convert  to  Cartesian 
coordinates  are 
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x  -  rcosQ 


y  -  rsinQ 


r-  yjx2*y2 


(65) 


Note  that  the  equation  for  G  above  will  introduce  an 
ambiguity. 


Case  1:  0  £  180*- y 


X  -  rcose  -  VFT7c„s[(I^l_)tan-'(^)] 
5-  rsin«-  V*2*  y2s‘n[(f5^)tan"^)] 


Case  2:  0>  180*-  y 


x 


y 


rcosQ  -  -V*2-*-  y2cos^~- j^tan'J^  180* 

rsin0  -  -Vx2+  y2sin^~-j^tan‘‘^  j  -  180°  jj 


(66) 


(67) 


(68) 


(69) 
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The  equations  for  the  rotation  are  considerably  more 

complicated  than  for  the  a  shift.  Not  only  does  each 
equation  have  complicated  terms ,  but  there  are  different 
equations  for  different  parts  of  the  plane.  Some  type  of 
'if /then'  programming  logic  would  be  required  to  implement 
them.  Also,  the  tan'1  term  in  each  equation  will  cause  a 
problem.  Mathematically,  two  solutions  are  possible,  but 
only  one  is  correct. 


4.3.  Combined  a  Shift  and  t,  Rotation 

The  next  logical  step  would  be  to  combine  an  a  shift  and 

a  £  rotation.  The  easiest  method  is  to  first  rotate  the 
imaginary  axis  to  some  angle  y  and  then  shift  the  entire 
plane  by  a.  Pigure  6  shows  these  transformations.  For  this 
section,  the  first  transformation  is  designated  by  a  bar  Cy) 
and  the  second  transformation  is  designated  by  a  hat  (?). 

The  equations  for  this  are  very  similar  to  the  ones  from 
Section  4.2.2  The  new  equations  are 

9my  £  -x  +  a  (70) 
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#.^T723m[(Ii^)tan”(j)] 


(72) 
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For  0  >  180°  -  y: 


-  )-  180*  )(+  a 

(73) 

LA  y  A  \> 

cJ  )\ 

(74) 


These  equations  are  not  much  more  complicated  than  the  t, 

rotation.  The  equations  for  the  £  rotation  and  the  a  shift 
can  be  combined  directly  with  no  additional  effort.  The 
next  section  will  show  how  the  transformations  affect  a 
vertical  line  at  -7. 

4.4.  Effects  of  the  Transformation  on  3  Vertical  Line 

He  will  begin  with  a  vertical  line  at  x  =  -7  and  first 
perform  just  a  t,  rotation,  and  then  both  a  t,  rotation  and  an 
a  shift.  Figure  7  shows  the  vertical  line  at  -7  and  also 
the  new  y  axis  after  a  t;  rotation  of  45°.  The  hatched  line 
shows  the  boundary  of  the  new  region  that  will  be  considered 
stable . 


30 


Figure  7  Line  before  coordinate  transformation 

Figure  8  shows  the  line  after  a  t  rotation  of  45°.  Note 

how  it  begins  at  x  =  -7  and  crosses  the  new  x  =  0  line  just 
below  y  =  10.  From  Figure  7  the  distance  to  the  point  where 
the  vertical  line  crosses  the  45®  line  is  9.89.  Since  r  is 
held  constant  over  the  transformation  this  is  exactly  what 
was  expected. 

Figure  9  shows  the  line  after  a  £  rotation  of  45®  and  an 
a  shift  of  5.  It  is  simply  Figure  8  shifted  by  5. 
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-7 


0 


2 


Figure  8  Line  after  fc;  rotation 


Figure  9  Line  after  combined  transformation 
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This  chapter  derived  the  equations  for  transforming  the 
complex  plane.  The  next  chapter  will  put  the  a  shift  into 
state  space  equations. 
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5.  State  Space  Equations  for  the  a  shift 

In  order  to  be  used  in  the  H2  and  H .  optimization 

equations  from  Chapter  2,  the  equations  derived  in  the 
previous  chapter  must  be  converted  from  Cartesian 
coordinates  to  state  space  equations.  The  easiest  way  to  do 
this  is  to  first  convert  them  to  the  Laplace  Domain  s-plane 
and  then  to  state  space.  For  the  a  shift  this  is  relatively 
easy.  The  y  (imaginary)  axis  is  moved  to  the  left  by  a  and 
the  x  (real)  axis  is  not  changed.  Looking  at  this  in  terms 
of  the  s-plane, 

s'-s  +  a  (75) 

or 

s - s ' -  a  (76) 

where  s'  denotes  the  transformed  s-plane. 

The  transfer  function  matrix  equation  relating  the  state 
space  equations  to  the  Laplace  domain  is 

G(s)-C(sI-AylB  <77> 


34 


Substituting  s»s'-a  into  this  equation  yields 


G(s')-  C((s'-a)/-  A)'1  B 


(78) 


•C(s'I-*I- AY'B 


(79) 


Let 


A' -al  +  A- A  +  al  (80) 

Substituting  this  into  (77)  yields 

G(s')-C(s'J-v*')'‘B  (81) 

This  is  the  standard  transfer  function  matrix  form 
again.  These  results  appear  rather  simple,  but  are  very 
powerful.  To  a  shift  the  state  space  equations,  all  that  is 
required  is  to  add  a /  to  the  A  matrix.  The  new  state  space 
equations  are 


x-  A’x+  Bu 


(82) 


yCx 


(83) 
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where 


A'  •  A  +  al  (84) 

As  described  in  Section  3.1,  these  equations  were  first 
derived  by  Anderson  and  Moore  in  1969  [3].  They  worked 
completely  in  the  time  domain  using  e*‘  as  a  weighting  factor 
on  the  performance  index.  Their  derivation  and  what  is 
presented  here  can  be  reconciled  by  the  Laplace  Transform 
identity 


e“F( 0  *  /(s-oO 


(85) 
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6.  Example  Problem  JL 


6.1.  Block  Diagrams  and  Generalised  State  Space  Squat icas. 

The  first  example  problem  uses  the  system  configuration 
shown  in  Figure  10.  The  outputs  of  interest  are  ei  and  e2. 
The  output  ei  weighted  by  W»  provides  control  over  the  shape 
of  the  sensitivity  function  and  the  output  e2  weighted  by  Hu 
provides  control  over  the  amount  of  control  power  the 
compensator  will  use. 


Figure  10  Block  diagram  for  Example  1  with  the  weighting 

factors 

The  first  step  is  to  put  the  system  into  standard  form. 
Standard  form  requires  state  space  equation  of  the  form 
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x  -  Ax*  Bdd*  Bau 


(86) 


e-C.x  +  D^d*  D^u 


(87) 


y~Cr+  Dydd  +  Dyttu 


(88) 


where  the  elements  of  the  above  equations  must  be  found  by 
manipulating  the  block  diagram  in  Figure  10. 

From  Pigure  10 


x - Apx  +  Bpu 


(89) 


yp-Cpx  +  Dpu-Cpx 


(90) 


y-d+y„-d+  Cpx 


(91) 


where  G  is  assumed  to  be  strictly  proper. 

The  state  space  equations  for  the  weighting  factor  W„ 


are: 


x#-  A,x,+  B,u ,  -  A,x,+  B,(d  +  C pX) 


(92) 


38 


(93) 


e,  -  C,x,+  Dtu,"  C,x, 

For  the  Doyle/Glover  parameterization  given  in  Chapter 
2,  D«d  must  be  equal  to  0.  To  insure  this,  we  require  that 
the  weighting  factor  Wt  be  strictly  proper,  i.e.  D»  =  0. 

The  reason  for  this  can  be  seen  from  Figure  10.  If  l/t  had 
a  Da  term,  then  the  disturbance  would  have  an  open  path  to 
the  output  ei  and  D«u  would  not  equal  0. 

In  this  example  we  will  restrict  the  weighting  l/„  to  be 
purely  a  gain,  i.e.  only  a  Du  term,  but  in  general  the  state 
space  equations  for  this  block  are 

+ (94) 

c2 “  Caxa+  Dau  (95) 

For  the  Doyle/Glover  parameterization,  D«uTD«u  must  be 
similar  to  I;  that  is,  it  must  be  full  rank.  In  practice, 
this  means  that  each  element  in  u  must  go  to  one  of  the 
controlled  outputs  directly.  To  insure  this,  the  l/K 
weighting  factor  must  have  a  Du  term.  Again  looking  at  the 
block  diagram  in  Figure  10,  the  reason  for  this  is  obvious. 
Since  G  and  Wt  are  strictly  proper,  ei  does  not  pick  up  u. 
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therefore  u  must  pass  through  W u  by  way  of  a  Du  term. 

Putting  this  in  matrix  form  and  dropping  the  Wu  state 
equations,  the  system  becomes 


(e‘)*(o  o)(*.)*(o  )rf*(o.)“ 


a 


y-(Cp  0)  :  U(7)d  +  (0)u 


(96) 


(97) 

(98) 


From  this  we  get  the  following  matrices  for  each  of  the 
variables  in  Equations  (86)  -  (88): 


*■■(*•) 

MS  e.O  »--(S) 

b“'(dJ 

C,-(CP  0)  £>*-(/) 

b„-(0) 

(99) 
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6.2.  Numerical  Values  for  the  System  Used 
The  test  plant  used  for  this  example  was 


G(s) 


2(s  +  2)(s  +  4) 
(s+  l)(s  +  3)(s  +  5) 


(100) 


This  transfer  function  was  chosen  to  show  the  results  of 
shifting  the  region  considered  to  be  stable  across  several 
open-loop  poles  and  zeros.  For  example,  an  a  shift  of  1.5 
would  require  a  compensator  designed  for  an  unstable  system 
in  the  s'  plane  and  an  a  shift  of  2.5  would  require  a 
compensator  designed  for  an  unstable  and  nonminimum  phase 
system  in  the  s'  plane. 

Converting  this  to  state  space  equations  gives 


C,-(  2  -4  -3.6148)  Dp-m 


(101) 


The  W ,  weighting  function  used  for  this  example  was: 
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1000(s+  1) 

(s  +  .001  )(s  +  1000) 


(102) 


V.- 


Figure  11  shows  a  singular  value  plot  of  W ,.  This  is 

chosen  to  have  a  shape  that  is  the  inverse  of  the  desired 
sensitivity.  As  can  be  seen  in  Figure  11,  the  l/,  weighting 
rolls  off  at  very  high  frequency  -  this  is  necessary  due  to 
the  strictly  proper  requirement.  The  desired  sensitivity 
does  not  look  like  the  inverse  of  W,  here,  but  since  these 
frequencies  are  well  beyond  the  bandwidth  of  the  system, 
this  will  not  affect  the  design.  Converting  this  to  state 
space  equations 

*-rr  v)  ■•-(;) 

C.-O000  1000)  £>,-(  0)  (103) 

The  Wu  used  was 

Dm  -  0.01  (104) 
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rad/ sec 


Figure  11  Singular  Value  plot  of  Wt 


Putting  all  this  together  into  the  standard  format 
yields : 
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(0  0  0  1000  1000\ 
Vo  0  0  0  0  J 


Cy  -  (2  -4  -3.6148  0  0) 


°--(o)  M.0°>) 


I)  Dya-(0)  (105) 

The  above  equations  are  in  standard  form  for  designing 
Hz  and  Hm  controllers.  To  close  the  loop  and  evaluate  the 
performance  of  the  closed-loop  system,  a  slightly  different 
form  is  required. 


6.3.  Closing  the  Loop 

The  block  diagram  configuration  used  to  close  the  loop 
is  shown  in  Figure  12.  This  is  still  standard  form.  To  get 
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this  configuration  the  block  diagram  shown  in  Figure  11  must 
be  modified  to  include  an  input  r.  This  modification  is 
shown  in  Figure  13. 


Figure  12  Block  diagram  for  Example  1  in  feedback  format 


Figure  13  Block  diagram  for  Example  1  with  command  input 


The  equations  for  this  figure  are 
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(106) 


x - Apx  + Bpu 

y  P-Cpx  (107) 

emyp  +  dmCpx  +  d  (108) 

y  -  r  +  e  -  r  +  Cpx  +  d  (109) 

Putting  these  equations  into  matrix  form  yields 

fd\  (HO) 

x-dpx  +  (0  0)[  J+  Bpu 

(:)-(">*(? 

/d')  (112) 

y -cpx  +  (/  /)(^J 

Putting  the  numbers  for  Example  1  into  the  above 
equations  yields 
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Cy  -  (2  -4  -3.6148) 


Dyd-(  1  1) 


Dyu  -  0  (113) 


6.4.  Mechanics  of  the  a  Shift 

To  accomplish  the  a  shift,  only  the  Ap  portion  of  the  A 

matrix  (see  Equation  99)  is  changed.  It  is  important  that 
a /  is  added  only  to  A„,  and  not  to  the  portion  of  the  A 
matrix  that  contains  the  weighting  functions.  The 
compensator  K  is  then  designed  using  standard  Hz  and  //. 
optimization  methods.  Once  the  compensator  is  designed  it 
must  be  reverse  shifted  back  to  the  standard  s-plane.  This 
is  accomplished  by  subtracting  a /  from  Ak,  where  Ak  is  the 
A  matrix  of  the  compensator. 
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6.5.  Results  for  Ho  Optimization 


The  next  step  is  to  perform  Hz  and  H m  compensator 

design  on  the  matrix  equations  derived  in  Section  6.2.  This 
will  yield  a  compensator  that  minimizes  the  2  or  »  norm  of 
the  weighted  closed-loop  system.  The  loop  is  then  closed 
around  the  compensator  using  the  block  diagram  and  the 
matrix  equations  derived  in  Section  6.3.  For  the  sake  of 
brevity,  only  selected  a  shifts  for  the  Hz  optimization 
design  method  are  presented  here.  Table  2  at  the  end  of 
Section  6.5  compares  the  closed-loop  poles  and  transmission 
zeros  for  a  =  0,  1.5  and  2.5.  The  shifted  plant  for  a  =  0 
is  stable  and  minimum  phase,  for  a  =  1.5  is  unstable  and  for 
a  =  2.5  is  unstable  and  nonminimum  phase. 

6.5.1.  a  -  0  (Original  Plant). 

This  section  shows  how  the  closed-loop  system  performs 
with  no  a  shift;  i.e.  standard  H2  optimization.  Figure  14 
shows  the  sensitivity  for  ue  system.  It  starts  at  -60  db, 
breaks  upward  at  w  =  0.001  and  levels  off  at  0  db.  This  is 
a  very  good  sensitivity  and  is  exactly  what  was  asked  for  by 
the  weighting  function  W». 
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Figure  14  Sensitivity  for  Example  1  with  no  a  shift 

Figure  15  shows  the  amount  of  control  energy  required. 
Note  the  dip  around  1  rad/sec  *nd  the  peak  just  below  20  db. 


rad/sec 

Figure  15  Control  usage  with  no  a  shift 
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Figure  16  shows  the  time  response  of  the  system  due  to  a 
unit  step  command  input  (e  to  r).  Mote  the  excellent 
command  following  and  the  relatively  slow  rise  time. 


Figure  16  Response  to  a  unit  step  command  input  with  no  a 

shift 

As  can  be  seen  from  the  above  figures,  this  was  a  very 
simple  case.  The  plant  was  already  stable  and  minimum 
phase.  The  sensitivity  has  a  good  shape  and  the  time 
response  tracked  to  the  desired  final  value  of  1  perfectly. 

6.5.2.  a  -  1 .5  (Unstable) . 

This  case  is  a  little  more  difficult.  Shifting  all  of 
the  poles  left  of  -1.5  will  require  "stabilising”  the  pole 
at  -1  from  the  plant.  Figure  17  shows  the  sensitivity. 
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Note  how  much  the  low  frequency  end  has  changed.  For  the 
unshifted  case  it  started  at  -60  db  (as  it  should  because  of 
the  gain  of  1000  in  N»),  but  for  the  shift  of  1.5  it  starts 
at  -12  db. 


Figure  17  Sensitivity  for  Example  1  with  1.5  a  shift 

Figure  18  shows  the  control  power  required.  Note  the 
slight  drop  at  low  frequencies  and  the  increase  around  1 
rad/sec. 

Figure  19  shows  the  time  response  to  a  unit  step  command 
step  input  (e  to  d).  The  system  responds  much  faster,  but 
the  final  value  is  significantly  degraded  to  0.7  rather  than 
1. 
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Figure  lb  Control  usage  with  1.5  a  shift 


Figure  19  Response  to  a  unit  step  command  input  with  1.5a 

shift 


6.5.3.  a ■  2.5  (Unstable  and  Nonminimum  Phase). 

Shifting  the  region  considered  to  be  stable  left  of  -2.5 
makes  the  problem  considerably  more  difficult.  Now  not  only 
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must  the  compensator  move  the  pole  at  -1,  but  also  must 
avoid  having  the  sero  at  -2  pull  a  pole  right  of  -2.5. 
Figure  20  shows  the  sensitivity.  Again  note  the  starting 
point  of  -12  db.  As  can  be  seen,  increasing  the  shift  from 
1.5  to  2.5  does  not  affect  the  sensitivity  very  much. 


rad/sec 

Figure  20  Sensitivity  for  Example  1  with  2.5  a  shift 

Figure  21  shows  the  control  power  required.  Note  the 
rather  large  increase  around  1  rad/sec. 

Figure  22  shows  the  time  response  to  a  unit  step  command 
input.  Note  the  greatly  increased  speed  of  response,  the 
final  value  of  1.2,  and  for  the  first  time  an  overshoot. 
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Figure  21  Control  usage  with  2.5a  shift 


Figure  22  Time  response  to  a  unit  step  command  input  for 
Example  1  with  a  2.5  a  shift 

6.5.4.  Conclusions  for  H*  Optimisation. 

Table  1  gives  a  suninary  of  the  results  shown  above,  as 
well  as  for  other  a  shifts.  Note  how  the  2  norm,  the 
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control  usage,  the  control  rate,  and  the  steady  state  error 
increase  as  a  is  increased.  Also  note  how  the  speed 
increases  up  to  a  certain  point  and  then  slows  down  again. 
Obviously,  only  the  smaller  a  shifts  are  reasonable  for  any 
kind  of  design. 

Table  1  Summary  of  H2  optimization  results 

a  shift 


0 

1.5 

2.5 

3.5 

5.5 

2  norm 

246.7 

249.4 

283.5 

501.5 

11514 

Control  usage 

5.8 

6.0 

7.6 

22 

610 

Control  rate 

3000 

3000 

4000 

1100 

61000 

Rejection  error 

0% 

30% 

30% 

80% 

200% 

Rejection  time 

2.5 

1 

0.5 

1 

1 

Tracking  error 

0% 

30% 

30% 

80% 

200% 

Tracking  time 

2.5 

1 

0.5 

1 

1 

Table  2  shows  the  closed-loop  poles  and  zeros.  As  can 
be  seen,  the  a  shift  does  exactly  what  it  is  designed  to  do; 
specifically,  it  moves  all  of  the  closed-loop  poles  left  of 
-a. 
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Table  2  Closed-loop  poles  and  zeros  for  Example  1 


a 

=  0 

a  = 

1.5 

a  = 

2.5 

poles 

zeros 

poles 

zeros 

poles 

zeros 

-978.9 

-11.78 

-980.4 

-24.16 

-981.4 

-380.2 

-204.3 

-5 

-205.8 

-5 

-206.8 

-5 

-5 

-5 

-5 

-5 

-5 

-5 

-4 

-3 

-4 

-3 

-4 

-3 

-3 

-3 

-3 

-3 

-4 

-2.662 

-2 

-1 

-2.5 

-1 

-3.5 

-1 

-1 

-1 

-2 

-3 

-1 

-2 

-3 

Another  interesting  result  is  the  change  in  the  2  norm. 
For  the  unshifted  case  the  2  norm  is  246.7.  For  an  a  shift 
of  1.5  the  2  norm  has  increased  to  249.4  and  for  an  a  shift 
of  2.5  the  2  norm  has  increased  to  283.5.  Note  that  there 
is  relatively  little  change  as  the  shift  crosses  a  pole,  but 
a  much  larger  increase  as  it  crosses  a  zero.  This  is  to  be 
expected.  Crossing  a  pole  is  the  same  as  putting  one  more 
unstable  pole  in  the  plant.  This  is  not  a  major  problem  for 
a  minimum  phase  system.  It  is  relatively  easy  to  move  poles 
around.  Crossing  a  zero,  however,  is  more  of  a  problem. 
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This  is  the  same  as  having  one  more  nonminimum  phase  zero  in 
the  plant.  The  poles  move  to  the  zeros  as  gain  is 
increased,  so  having  a  pole  and  a  zero  in  the  right  half 
plane  is  a  very  tricky  problem.  The  compensator  must  send  a 
pole  into  the  right  half  plane  to  split  the  unstable  plant 
pole  off  the  real  axis  and  pull  them  both  into  the  left  half 
plane  without  allowing  the  zero  to  draw  a  pole  to  the  right 
half  plane. 

The  shift  of  the  closed-loop  poles  shows  up  very  clearly 
in  the  time  response  plots.  As  a  is  increased  the  time 
response  get  much  faster.  One  problem  is  the  degradation  of 
the  final  value  that  the  time  response  approaches.  For  good 
command  following  it  should  go  to  1.  One  way  to  improve 
this  would  be  to  add  a  bank  of  integrators  to  the  system. 
Another  possible  solution  would  be  to  experiment  with  the  W» 
weighting  factor  or  to  decrease  the  Hu  parameter  to  allow 
the  system  to  use  more  control  power. 

6.6.  RgguLfcg  tQ*  H-  £>ptimi*atign 

The  shape  of  the  plots  and  the  effects  of  the  a  shift 

are  very  similar  to  H2  optimization  for  Hm  optimization. 

The  magnitude  does  change  in  the  plots,  but  trends  remain 
the  same.  The  results  are  summarized  in  Table  3. 
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Example  1  showed  the  results  of  shifting  across  a  pole 
and/or  zero,  and  the  trends  and  limits  as  a  is  increased. 
It  also  presented  the  mechanics  of  the  method.  The  next 
example  will  test  this  method  on  a  real  design  problem. 

Table  3  Summary  of  Hm  optimization  results 

a  shift 


0 

1.5 

2.5 

3.5 

5.5 

Gamma 

0.855 

0.856 

6.009 

13.44 

268 

Control  usage 

75 

75 

550 

1200 

25000 

Control  rate 

75000 

75000 

5.5*105 

1 . 2*10® 

2.5*107 

Rejection  error 

0% 

25% 

100% 

120% 

200% 

Rejection  time 

2.5 

1 

0.5 

0.5 

1 

Tracking  error 

0% 

25% 

100% 

120% 

200% 

Tracking  time 

2.5 

1 

0.5 

0.5 

1 
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7.  Example  Problem  2. 


7.1.  Problem  Statement 

The  previous  example  showed  the  mechanics  of  how  the  a 

shift  works.  This  example  applies  it  to  a  problem  taken 
from  [10:13-1].  The  problem  is  to  design  a  lateral  attitude 
control  system  for  a  drone  aircraft.  The  state  space 
equations  for  this  example  are 


x - Ax  + Bu 


(114) 


y - Cx + Du 


(115) 


where : 


A 


/- 0.08527 
-46.86 
-  0.4248 
0 

0° 


-0.0001423 

-0.9994 

0.04141 

-  2.757 

0.3896 

0 

-  0.0622 4 

-0.6714 

0 

1 

0.0523 

0 

0 

0 

0 

0 

0 

0 

0  0.1862  \ 

-124.3  128.6  \ 

-8.792  -  20.46 

0  0 

-20  0 

0  -20  / 


(116) 
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f°0 

B- 

1  0 

0 

!  0 

0 

[  20 

0 

(117) 

\  o 

20/ 

c-f1 

0  0 

0  0 

(118) 

C  lo 

0  0 

1  0 

o) 

n 

m(° 

°) 

(119) 

Is 

Vo 

oj 

x-[P  i 

V 

♦  6. 

MT 

(120) 

u- 

IK 

Orcf 

(121) 

The  states  x  are  sideslip  angle,  roll  rate,  yaw  rate, 
roll  angle,  elevon  surface  deflection,  and  rudder  surface 
deflection.  The  inputs  u  are  elevon  servo  command  and 
rudder  servo  command.  The  open- loop  plant  has  one  zero  at 
-158.15.  The  open-loop  poles  are  shown  in  Table  4. 
Obviously,  the  compensator  must  stabilize  the  unstable  dutch 
roll  mode.  The  following  figures  and  tables  show  that  Hz 
optimization  simply  pulls  that  complex  pair  to  its  mirror 
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image  on  the  left  side  of  the  imaginary  axis  and  then  adds 
complex  zeros  very  close  by  to  nearly  cancel  the  effects  of 
the  lightly  damped  complex  poles. 


Table  4  Open- 

loop  poles 

Dutch  Roll  (unstable) 

0.1884  a  1.0511  i 

Spiral  Mode 

-0.0360 

Roll  Convergence 

-3.2503 

Elevon  Actuator 

-20.0 

Rudder  Actuator 

o 

• 

o 

<N 

1 

7.2.  H ,  Optimization  Setup  and  Equations 

Figure  23  shows  the  block  diagram  for  designing  the 
compensator.  Mote  the  additional  output  ( e2 )  which  is 
simply  the  output  from  the  plant.  The  output  e2  is  required 
for  this  problem  because  we  want  to  control  the 
complimentary  sensitivity  as  well  as  the  sensitivity.  The 
complimentary  sensitivity  is  the  transfer  function  from  d  to 
e?.  The  other  two  outputs,  C)  and  e3,  are  the  same  ones  as 
used  in  the  previous  example.  e3  is  the  control  input  u 
weighted  by  an  identity  matrix  multiplied  by  some  constant 
p.  The  parameter  p  is  the  weighting  penalty  on  control 
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usage.  As  p  is  lowered  the  system  is  allowed  to  use  more 
control  power.  Por  this  problem  p  was  chosen  to  be  0.05. 
For  p  =  0.05  the  maximum  deflections  and  rates  of  the 
control  surfaces  (rudder  and  elevon)  are  within  the  range 
allowed  on  a  real  aircraft. 


y 


Figure  23  Block  diagram  for  design  on  Example  2 

e,  is  the  output  y  weighted  by  some  strictly  proper  transfer 

function  matrix  W W ,  is  chosen  so  that,  for  low 
frequencies,  the  sensitivity  will  be  its  mirror  image  across 
the  0  db  line.  In  Example  2,  W,  is  a  2x2  diagonal  transfer 
function  matrix.  The  1,1  and  2,2  elements  are  the  same  and 
equal  to 
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w 


1000($+ 1) 

($+  .001)(s  +  1000) 


i-  1.2 


(122) 


•u 


This  is  the  same  as  the  W,  from  Example  1.  In  the 
following  equations  the  subscript  p  refers  to  the  plant  (G) 
and  the  subscript  s  refers  to  the  I/a  weighting  function. 

The  equations  we  are  working  towards  are 


x  -  Ax*  Bdd+  Btu 

(123) 

e-  C,x+  D^d*  Dtau 

(124) 

y-C,  +  Dytd+Dy.u 

(125) 

Following  the  same  procedure  as  in  Example  1  the  final 
state  space  equations  become 


(126) 


(127) 
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(128) 


y-(  cp 


+  (/)d  +  (0)u 


The  above  setup  is  used  only  to  find  the  compensator  K. 
The  closed-loop  format  is  the  same  as  from  Example  1. 
Figure  24  shows  the  standard  HIMO  feedback  configuration. 
As  in  Example  1,  the  block  diagram  in  Figure  24  must  first 
be  modified  to  include  a  command  input  r.  The  weighting 
factors  and  W ,  must  also  be  removed.  They  are  included 
only  as  design  parameters  to  find  the  desired  compensator. 
Figure  25  shows  the  modified  block  diagram  used  to  find 
Figure  24. 


Figure  24  Block  uiagram  for  Example  2  in  feedback  format 
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Figure  25  Block  diagram  for  Example  2  with  command  input 


The  state  space  equations  for  this  for*  are 


x  m  A px  +  (0  0) 


(?)**’“ 


(129) 


CMcX?  XK> 


(130) 


y-Cpx  +  C/  /) 


(?) 


(131) 


These  equations  will  be  used  to  test  the  response  of  the 
system  to  both  a  step  disturbance  and  a  step  command  input. 
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7.3.  H2  Optimisation  with  ns.  a  Shill 


Por  p  =  0.05  and  no  a  shift  the  2  norm  is  12.12.  The 

sensitivity  is  shown  in  Figure  26.  This  sensitivity  is 
exactly  what  was  asked  for  by  the  weighting  function  W ,. 


rad/ sec 

Figure  26  Sensitivity  for  Example  2  with  no  a  shift 

Figure  27  shows  the  control  surface  deflection  in 
degrees  required  for  a  unit  step  $  (sideslip  angle)  command. 
The  maximum  deflection  is  slightly  more  than  -0.6  degrees 
and  the  maximum  rate  is  less  than  15  degrees  per  second. 
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Pigure  27  Control  power  required  for  a  unit  step  {S  command 

Figure  28  shows  the  control  surface  deflection  in 
degrees  required  for  a  unit  step  4>  (roll  angle)  command. 

The  maximum  deflection  is  less  than  0.4  degrees  and  the 
maximum  rate  is  less  than  20  degrees  per  second. 

Pigure  29  shows  how  the  system  responds  to  a  unit  step 
sideslip  angle  disturbance.  In  about  10  seconds  the 
disturbance  is  completely  rejected.  The  sideslip  returns  to 
0  fairly  quickly  with  minimal  oscillation. 
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Figure  30  shows  how  the  system  responds  to  a  unit  step 
roll  angle  disturbance.  Mote  the  time  scale  on  this  figure. 
The  oscillations  extend  beyond  20  seconds. 


time 


Figure  30  Rejecting  a  unit  step  $  disturbance 


Figure  31  gives  the  command  response  of  the  system  to  a 
unit  step  $  command.  Within  5  seconds  the  system  goes 
exactly  to  the  desired  value. 

Figure  32  gives  the  command  response  of  the  system  to  a 
unit  step  4>  command.  Note  the  time  scale  again.  The 
response  still  has  oscillations  at  20  seconds. 
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time 


to  a  unit 


step  ^  command  inpu 


T 


<P 


P 


20  30 

time 


to  a  unit  step  $  command  input 


Overall,  the  system  performance  is  not  too  bad.  It  will 
completely  reject  any  disturbance  and  tracks  command  inputs 
perfectly  in  the  steady  state.  The  control  deflections  and 
rates  are  well  within  the  capabilities  of  the  actuators. 

The  (3  disturbance  rejection  and  command  following  are  fairly 
quick  and  have  little  oscillation.  The  only  problem  is  in 
the  4>  disturbance  rejection  and  command  following.  The 
oscillations  do  not  damp  out  for  at  least  20  seconds.  Table 
5  shows  the  low  frequency  closed-loop  poles  and  zeros  for 
the  system.  There  are  an  additional  12  poles  and  8  zeros 
with  a  magnitude  greater  than  10.  They  are  not  shown  here 
because  they  have  very  little  effect  on  the  response. 

Figure  33  shows  the  location  of  the  poles  and  zeros  with 
a  magnitude  less  than  10  in  the  s-plane. 

The  open- loop  poles  at  0.18842  a  1.0511i  have  been  moved 

to  their  mirror  imag*  location  in  the  left  half  plane.  The 
zero  at  -0.057927  a  0.90067i  is  close  enough  to  make  the 
residue  from  the  poles  fairly  small,  but  the  poles  have  such 
poor  damping  that  they  cause  the  ringing  effect  in  Figures 
30  and  32. 
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Table  5  Closed-loop  poles  and  zeros 


poles 

zeros 

-0.035987 

0.18842  *  1.0511  i 

-0.77686  *  0.32181  i 

-0.035987 

-0.77689  *  0.32180  i 

-0.035987 

-0.18842  *  1.0511  i 

-0.057927  *  0.90067  i 

-3.2503 

-0.41344 

-0.43689 

-3.2503 

-3.2503 

■* — t— 
-3 


~-2 


o  zero 

•  double  zero 
x  pole 

*.  double  pole 


* 


4— 

-1 


* 


Figure  33  Poles  and  zeros  for  the  unshifted  case 
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The  next  section  will  show  how  an  a  shift  of  0.3  changes 


the  damping  on  those  poles,  and  therefore  the  time  response. 
As  will  be  seen,  the  cost  of  the  increased  damping  is 
decreased  performance  elsewhere. 

7.4.  H 2  Optimisation  with  AH  a  SM11  fit  0 ■  3 

The  mechanics  of  the  a  shift  are  the  same  as  for  Example 

1.  Only  the  Ap  portion  of  the  A  matrix  is  shifted,  and  Ax 
must  be  reverse  shifted.  An  a  shift  of  0.3  seemed  to  give 
the  best  trade-off  between  increased  damping  and  decreased 
performance.  This  moves  the  imaginary  axis  to  -0.3  and  the 
unstable  poles  at  0.18842  *  l.OSlli  are  moved  to  their 
mirror  image  across  the  new  imaginary  axis.  Their 
closed-loop  position  is  -0.78842  *  1.0511i.  Table  6  shows  a 
summary  of  the  results  of  an  a  shift  from  0  to  0.9.  Por  p  = 
0.05  and  an  a  shift  of  0.3  the  2  norm  ir  17.52.  The  2  norm 
for  the  unshifted  case  was  12.12.  As  should  be  expected, 
the  2  norm  has  increased.  The  compensator  is  required  to  do 
more  than  before,  and  it  will  not  be  able  to  get  to  the  same 
minimum  2  norm. 


73 


Table  6  Summary  of  H 2  optimization  results 


a  shift 


2  norm 
(3  dist. 

3  dist. 
<t>  dist. 
4>  dist. 
(3  track 
P  track 
4>  track 
<)>  track 


error  (%) 
time  (sec) 
error  (%) 
time  (sec) 
error  (%) 
time  (sec) 
error  (%) 
time  (sec) 


Figure  34  shows  the  sensitivity.  Note  the  changes  at 
low  frequency.  Before  we  had  both  channels  beginning  below 
-100  db  and  going  smoothly  to  zero  db  at  1  rad/sec.  For  the 
shifted  case  one  channel  begins  just  below  -20  db  and 
neither  of  them  are  exactly  what  we  wanted. 
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10-3  10°  103 

rad / sec 


Figure  34  Sensitivity  for  Example  2  with  0.3  a  shift 

Figure  35  shows  the  control  surface  deflections  required 
to  perform  a  unit  step  $  (sideslip  angle)  command  input. 

The  peak  deflections  and  rates  are  similar  to  the  unshifted 
case . 

Figure  36  shows  the  control  surface  deflections  required 
to  perform  a  unit  step  $  (roll  angle)  command  input.  Here 
the  peak  deflection  has  increased  to  0.8  degrees  (from  0.4) 
and  the  rate  has  increased  to  around  30  degrees  per  second 
(from  20).  This  is  still  reasonable,  but  it  is  a  lot 
greater  than  the  unshifted  case.  The  system  must  perform 
faster  so  more  control  power  is  required. 
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Figure  37  shows  how  the  system  responds  to  a  unit  step  $ 

disturbance.  This  is  where  the  penalty  for  the  a  shift 
begins  to  become  apparent.  The  system  responds  much  faster 
and  the  damping  is  increased,  but  it  does  not  completely 
reject  the  disturbance.  The  steady  state  value  is  about 
0.08  rather  than  the  desired  value  of  sero. 

Figure  38  is  the  $  disturbance  response.  This  is  a 

vast  improvement  over  the  unshifted  case.  The  oscillations 
damp  out  in  less  than  10  seconds  (compared  to  30)  and  the 
response  has  no  steady  state  error. 


time 

Figure  37  Rejecting  a  unit  step  $  disturbance 
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1 


time 

Figure  38  Rejecting  a  unit  step  $  disturbance 

Figure  39  shows  the  response  to  a  unit  step  (3  command 

input.  Again  the  speed  of  the  system  has  increased  over  the 
unshifted  case,  but  the  steady  state  has  an  error  in  it. 

The  final  value  is  about  0.92  rather  than  1. 

Figure  40  is  the  $  response.  Here  there  is  a  major 

improvement  over  the  unshifted  case.  The  oscillations  damp 
out  in  about  8  seconds  (compared  to  30)  end  the  system  has 
no  steady  state  error. 
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System  response  to  a  unit  step  4>  command  input 
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Table  7  shows  the  low  frequency  poles.  Again  the  system 


has  an  additional  12  poles  and  8  zeros  beyond  -10. 


Table  7  Closed-loop  poles  and  zeros 


poles 

zeros 

-0.56401 

0.18842  *  1.0511  i 

-1.0769  *  0.32181  i 

-0.035987 

-1.0769  *  0.32180  i 

-0.36400  *  0.90850  i 

-0.78842  *  1.0511  i 

-0.52713  *  0.13156  i 

-3.2503 

-0.73613 

-3.2503 

-3.2503 

In  Figure  41  the  positions  of  the  poles  and  zeros  are 
shown  graphically.  Note  how  the  poorly  damped  poles  in  the 
unshifted  case  have  moved  much  farther  into  the  LHP. 

The  open-loop  poles  at  0.18842  *  1.05111  have  been  moved 

to  their  mirror  image  location  across  the  shifted  imaginary 
axis  at  -0.3  to  -0.78842  *  1.05111.  This  greatly  increases 
the  damping  of  the  system. 
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Figure  41  Poles  and  seros  for  an  a  shift  of  0.3 

One  other  interesting  point  that  came  from  Example  2  was 
the  lack  of  dependance  on  p.  Even  when  p  was  lowered  to 
0.0001,  giving  the  system  virtually  unlimited  control  power, 
the  steady  state  error  did  not  improve  much  for  the  shifted 
case.  The  degradation  in  performance  after  the  a  shift  was 
not  due  to  a  lack  of  control  pwwer  available,  but  rather  due 
to  the  mathematics  of  the  controller  design. 
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8.  Perturbing  the  Plant 


Example  2  shows  how  H2  compensator  design  moves  unstable 

complex  poles  to  their  mirror  image  location  across  the 
imaginary  axis,  and  then  places  seros  near  them  to  limit  the 
effect  they  will  have.  This  process  works  well  if  the 
system  model  is  accurate  and  does  not  change.  Systems  do, 
however,  change  with  time  and  conditions,  and  no  model  is 
perfect.  If  the  poles  are  not  where  the  compensator  expects 
them  to  be,  the  xeros  will  not  cancel  them,  and  the  effects 
of  these  pole  will  increase.  This  chapter  will  look  at  the 
problem  from  Example  2  with  two  different  perturbations  on 
the  A  matrix.  These  perturbations  may  come  from  flying  in  a 
different  flight  condition  (new  altitude,  different  weight 
due  to  burned  fuel,  etc)  or  from  modeling  errors.  In  the 
perturbed  cases,  the  compensator  was  designed  for  the 
original  plant  and  then  the  loop  was  closed  around  the 
perturbed  plant.  The  original  A  matrix  is 


/ -0.085 27 

-0.0001423 

-0.9994 

0.04141 

0 

0.1»62 

-46.86 

-2.757 

0.3896 

0 

-  124.3 

128.6 

-0.4248 

-0.06224 

-0.6714 

0 

-  8.792 

-20.46 

0 

1 

0.0523 

0 

0 

0 

0 

0 

0 

0 

-20 

0 

\  o 

0 

0 

0 

0 

-20 

(132) 
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Table  6  shows  the  eigenvalues  of  the  original  open- loop 
plant . 


Table  8  Open-loop  poles 


Dutch  Roll  (unstable) 

0.1884  *  1.0511  i 

Spiral  Mode 

-0.0360 

Roll  Convergence 

-3.2503 

Elevon  Actuator 

i 

N> 

O 

• 

o 

Rudder  Actuator 

© 

• 

o 

ts 

1 

8.1.  Perturbation  Humber  1 

The  first  perturbed  A  matrix  was 


/ 


A  - 


0.082610 

-46.86 

0.020024 

0 

0 

O 


-0.00021019 

-  0.9994 

0.04141 

-4.0>  ■'«* 

0.3896 

0 

-  0.06224 

-0.6714 

0 

1 

0.0523 

0 

0 

0 

0 

0 

0 

0 

0  0.1862  \ 

-  124.3  128.6  1 

-8.792  -20.46 

0  O 

-20  0 

O  -20  / 


(133) 


This 


changed  the  open- loop  poles  to  the  ones  listed  in 


Table  9. 
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Table  9  Open- loop  poles 


Dutch  Roll  (unstable) 

0.07223  a  1.0727  i 

Spiral  Mode 

-0.024969 

Roll  Convergence 

-4.2800 

Elevon  Actuator 

i 

M 

o 

• 

o 

Rudder  Actuator 

-20.0 

8.1.1.  Perturbation  Number  i  with  no  a  Shift. 

The  sensitivity  changed  significantly  from  the 
unperturbed  case,  but  it  is  still  acceptable.  The  two 
channels  are  no  longer  together  and  the  upper  singular  value 
increased  to  above  -100  db  for  low  frequencies.  The  2  norm 
has  increased  to  13.67  from  12.12  in  the  unperturbed, 
unshifted  case.  The  control  usage  is  very  similar  to  the 
unperturbed  case,  but  with  a  little  more  oscillation.  The 
roll  angle  (4>)  disturbance  rejection  and  command  following 
are  also  similar  to  the  unperturbed  case.  The  big  change  is 
in  the  sideslip  O)  disturbance  rejection  and  command 
following.  They  both  go  to  an  initial  deflection  of  3.5 
(unit  step  input)  and  oscillate  for  almost  40  seconds.  See 
Figures  42-48. 
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Pigure  45  Rejecting  a  unit  step  0  disturbance 
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System  response  to  a  unit  step  JJ  command  input 
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Figure  48  System  Response  to  a  unit  step  $  command  input 

8.1.2.  Perturbation  Humber  1  with  Q.3  a  Shift, 

Here  the  sensitivity  has  degraded  to  about  -5  db  at  low 
frequencies.  The  2  norm  is  almost  the  same  as  the 
unperturbed  case  with  0.3  a  shift.  It  has  increased  from 
17.52  to  17.53.  Again,  control  usage  is  very  similar  to  the 
unperturbed  case.  As  in  the  perturbed,  unshifted  case  the  $ 
disturbance  rejection  and  command  following  are  very  good. 
The  big  problem  is  in  the  (3  steady  state  error.  The  huge 
deflections  from  the  unshifted  case  have  been  largely 
removed,  but  the  steady  state  values  for  the  unit  step  ^ 


disturbance  are  around  0.2  and  -0.4.  The  steady  state 
values  for  a  unit  step  ^  command  are  0.8  and  0.4.  The 
response  reaches  the  final  values  for  both  $  disturbance 
rejection  and  0  command  following  much  faster  (around  10 
seconds)  but  the  final  values  are  not  what  we  want  (see 
Figures  49-55) 


rad/sec 


Figure  49  Sensitivity  of  perturbation  1  with  0.3  shift 
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8.2.  Perturbation  Humber  2 

The  results  for  the  first  perturbation  are  questionable 
as  to  whether  the  a  shift  is  an  improvement.  The  system 
does  damp  out  much  faster,  but  the  steady  state  values  are 
very  poor  in  some  cases.  Consider,  however,  what  happens 
when  the  plant  is  perturbed  to  a  different  value.  The 
following  is  the  second  perturbed  plant 


-0.082610 

-0.00021019 

-  0.9994 

0.04141 

0 

0.1862\ 

-  46.86 

-  5.781  S 

0.3896 

0 

-  124.3 

128.6 

-  0.020024 

-  0.06224 

-0.671 4 

0 

-  8.792 

-20.46 

0 

1 

0.0523 

0 

0 

0 

0 

0 

0 

0 

-20 

0 

0 

0 

0 

0 

0 

-20  / 
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The  eigenvalues  for  this  plant  are  listed  in  Table  10. 


Table  10  Open- loop  poles 


Dutch  Roll  (unstable) 

0.006979  *  0.91767  i 

Spiral  Mode 

-0.024754 

Roll  Convergence 

-5.9205 

Elevon  Actuator 

-20.0 

Rudder  Actuator 

-20.0 
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This  is  not  much  of  a  change  from  the  first 
perturbation,  but  the  closed-loop  system  with  no  a  shift  is 
now  unstable.  Any  input  or  disturbance  results  in  an 
unstable  output.  Using  an  a  shift  of  0.3,  however,  results 
in  a  stable  closed-loop  system.  The  results  are  very 
similar  to  the  shifted  case  of  perturbation  number  1.  They 
may  have  poor  disturbance  rejection  and  tracking  in  some 
channels,  but  they  are  definitely  better  than  the  unstable 
results  for  the  unshifted  case.  Table  11  shows  a  summary  the 
effects  of  perturbing  the  A  matrix. 
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Table  11 


Perturl 

a  shift  0 


2  norm 

13.67 

3  diet,  error  (%) 

0 

3  diet,  time  (sec) 

50 

4>dist.  error  (%) 

0 

4>  dist.  time  (sec) 

25 

3  track  error  (%) 

0 

3  track  time  (sec) 

50 

<t>  track  error  (%) 

0 

4>  track  time  (sec) 

25 

j  the  A  matrix 

ion  1  Perturbation  2 

0.3  0  0.3 


9.  Results  and  Conclusions 


Overall,  the  a  shift  did  exactly  what  it  was  designed  to 

do.  It  shifted  all  the  closed-loop  poles  left  of  a  vertical 
line  some  distance  a  from  the  imaginary  axis.  It  also 
increased  the  damping  of  complex  poles  very  near  the 
imaginary  axis.  There  were,  however,  significant  reductions 
in  the  steady  state  performance. 

9.1.  Summarising  the  Example  Problems 

Example  1  showed  how  the  system  responded  as  the  a  shift 

crossed  each  pole  or  xero.  Shifting  across  a  pole  (making 
the  system  "unstable")  or  xero  (making  the  system 
"nonminimum  phase")  caused  a  major  decrease  in  the 
performance.  The  low  frequency  sensitivity  increased,  the 
control  usage  increased,  and  the  steady  state  error  for 
disturbance  rejection  and  command  following  increased.  This 
was  reflected  in  drastic  increases  in  both  the  2  and  «  norm. 
Damping  was  not  a  factor  in  Example  1,  but  it  did  show  the 
desired  increase  in  the  speed  of  response. 

Example  2  was  a  more  interesting  test  of  the  method.  In 
this  example  there  were  very  lightly  damped  poles  in  the 
open-loop  system.  The  a  shift  did  considerably  increase  the 
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closed-loop  damping  of  these  poles.  The  open-loop  plant  was 
unstable,  with  a  complex  pair  just  right  of  the  imaginary 
axis.  Standard  H2  optimization  moved  these  poles  to 
roughly  their  mirror  image  locations  in  the  left  half  plane. 
This  stabilized  the  system,  but  these  poles  had  very  poor 
damping.  The  time  response  had  oscillations  for  more  than 
20  seconds. 

H2  optimization  with  an  a  shift  moved  the  open-loop 
unstable  poles  to  their  mirror  image  locations  across  the 
shifted  imaginary  axis.  In  effect  this  gave  double  the 
shift.  The  a  shift  dramatically  increased  the  damping  of 
the  closed-loop  system.  The  time  response  settled  out  in 
less  than  10  seconds.  The  cost  for  this  increased  damping 
was  mostly  in  the  steady  state  error.  There  was 
approximately  an  eight  percent  error  in  both  disturbance 
rejection  and  command  following  for  one  of  two  inputs.  The 
other  input  still  had  no  steady  state  error.  An  additional 
cost  was  a  minor  increase  in  the  control  power  required. 

Chapter  8  showed  how  the  a  shift  increased  the 
robustness  of  the  closed-loop  system.  The  A  matrix  of  the 
plant  was  perturbed  to  produce  two  new  A  matrices, 
representative  of  two  off-nominal  fight  conditions.  The 
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compensator  designed  using  an  cl  shift  was  able  to  stabilize 
larger  variations  in  the  A  matrix  than  standard  Hz 
optimization. 

9.2.  Conclusions 

The  a  shift  is  basically  another  design  variable,  like  p 

and  W It  is  most  useful  for  systems  that  have 
closed-loop  poles  near  the  imaginary  axis.  It  will  increase 
the  speed  of  response  of  the  closed-loop  system,  and  for 
lightly  damped  poles  near  the  imaginary  axis  it  will 
increase  the  damping  ratio  significantly.  The  a  required  is 
dependant  on  the  particular  system.  Due  to  the  tendency  of 
Hz  and  //_  optimization  to  move  unstable  poles  to  their 
mirror  image,  the  benefit  from  an  a  shift  is  double  what 
would  be  expected.  The  penalty  for  using  an  a  shift  can  be 
significant.  Moving  all  the  poles  away  from  the  origin 
removes  any  integral  action  that  may  have  been  in  the 
system.  The  steady  state  error  can  increase  dramatically, 
especially  as  a  get  large. 

9.2.1.  Areas  ion  Further  Study. 

The  first  thing  to  do  would  be  to  reproduce  the  results 
from  Example  2  using  optimization.  Most  likely,  they 
would  be  similar  to  the  H2  case,  buc  may  be  worth 
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considering . 

Next,  it  would  be  very  beneficial  to  get  the 
sensitivity,  disturbance  rejection,  and  command  following  of 
the  unshifted  case  to  carry  over  to  the  shifted  case.  There 
are  several  possible  methods  for  doing  this.  The  easiest 
things  to  try  would  be  to  adjust  the  weighting  factors  W , 
and  l/„.  Increasing  the  gain  on  W,  and  extending  the 
bandwidth  to  around  10  rad/sec  may  help.  For  V/H  it  might 
be  worth  looking  at  a  diagonal  weighting,  penalizing  the 
control  usage  differently  in  each  channel.  If  neither  of 
these  work,  one  could  consider  adding  an  additional  outer 
loop  with  integral  action. 

The  most  important  follow  on  work  would  be  to  derive  the 
state  space  equations  for  the  t;  rotation.  This  would  move 
the  lightly  damped  poles  away  from  the  imaginary  axis 
without  moving  the  poles  near  the  origin  that  cause  the  loss 
of  integral  action.  It  is  unlikely  that  the  t,  rotation  will 
work  out  as  simply  as  the  a  shift,  but  it  would  be  worth  the 
effort . 

Finally,  a  combined  shift  of  the  form  shown  in  Figure  56 
could  be  derived.  The  combined  transformation  shown  in 
Chapter  4  does  not  guarantee  a  specific  damping  ratio.  The 
new  y  axis  is  no  longer  rotated  about  the  origin,  but  about 
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some  point  a  from  the  origin.  The  shift  in  Figure  56  would 
increase  the  complexity  of  the  transformation,  however, 
because  only  a  portion  of  the  rotated  axis  is  shifted,  and 
the  shift  is  not  constant  for  all  values  of  y. 


Figure  56  Combined  s-plane  transformation 
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